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"64 ! Abstract 
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\ We reanalyze the conditions under which we have a primordial thermal population of axions. We 

extend a previous study and take into account processes involving gluons and quarks. We conclude 

a : 

that if the Peccei-Quinn scale fulfills Fa < 1.2 x 10^^ GeV there is thermal axion production. In 

tJ" ' this case, a period in the early universe exists where axions would interact with the QCD plasma 
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and we point out that non-thermal axions produced before the end of this period would thermalize. 
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I. INTRODUCTION 



The Peccei-Quinn solution [||] to the strong CP-problem introduces a new U{1)pq chiral 
symmetry, the spontaneous breaking of which leads to a new spinless particle, the axion 
0. If it exists, the axion has quite well defined properties: its mass and couplings are 
inversely proportional to Fa, the energy scale of the U{1)pq breaking. We will be particularly 
concerned with the coupling of the axion field a to gluons 

^gga = — ^ G^^" a (1) 

r a oTT 

where = g'^/^iir) with gs the gauge coupling of color SU{3)c. In (0), 

is the dual of the gluon field Gj^j,, and there is a sum over the color index b. 

Much theoretical study and experimental effort have been devoted to this particle (for 
recent reviews see [^). The laboratory, astrophysical, and cosmo logical constraints imply 
that Fa must be a high energy scale, at least as large as [§] 

Fa>6x 10® GeV (3) 

which in turn implies 

rria < 0.01 eV (4) 

Thus, the axion has to be a very light particle and very weakly coupled. 

One of the most attractive features of the axion is that it may be one of the constituents 
of the dark matter. Indeed, axions are copiously produced in the early universe. These 
primordial axions are generated mainly non-thermally. We know three non-thermal sources 
at work in the evolution of the universe: the "misalignment mechanism" [Q, string decay 
^ and domain wall decay [0]. When the Peccei-Quinn symmetry breaks down at a tem- 
perature T ~ Fa, the vacuum angle has no preferred value; but when the temperature 
reaches T ^ Aqcd, the vacuum angle sets where the vacuum energy is minimized, which 
is the CP-conserving value. Relaxation from one angle to the other produces a coherent, 
cold condensate of axions. This is the "misalignment mechanism" . On the other hand, the 
U{l)pQ breaking at T ~ originates axion strings. If inflation takes place just after that, 
these strings dilute away. However, if this is not the case, we have a second mechanism of 
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production of axions when strings decay. Also, if the axion field is not homogeneized by 
inflation, there are models where axion domain walls appear soon after T ~ ^qcd- The 
decay of these walls into axions is a third mechanism of non-thermal primordial axions. 

Thermal production of axions was studied by Turner |]^. His main interest was to study 
the possibility that thermal axions have a density greater than the corresponding to non- 
thermal axions, and then to analyze the potential detection of photons from axion decay. 
Turner considered the Primakoff process, a + g 7 + g, where g is a light quark, and 
photoproduction, a + Q ^ 7 + Q, where Q is a heavy quark, and showed that when 

Fa < 10^ GeV (5) 

axions were once in thermal equilibrium. Evaluating the axion density he found that for 

F„ < 4 X 10^ GeV (6) 

the thermal population is greater than the non-thermal. 

In [H, it was mentioned that other axion processes contribute to axion thermalization 
and consequently would increase the value in (||). However, these other processes were not 
evaluated since the range includes (H) and this was enough for the purposes of 0. We 
note that the experimental bound in @ excludes the range in (P) and leaves a small interval 
of the condition (j^). 

In this paper we extend the work of [§] and investigate the contribution to axion ther- 
malization of other processes than Primakoff and photoproduction. As we will show, the 
processes we discuss lead to a much larger range for Fa than (|^). 

The reactions we consider all involve the gga vertex appearing in (|I[), 

1) a + g^q + q 

2) a + q ^ g + q and a + q^ g + q 
S)a + g^g + g 

We show the corresponding Feynman diagrams in Figs. 0, |], and |^. 

The coupling of the axion to gluons (|l|) is an essential ingredient of the effective theory 
once the Peccei-Quinn symmetry is broken, since it reproduces the chiral anomaly of the 
theory. In this sense our results are model independent in contrast with the result @ arising 
from Primakoff and photoproduction processes that involve the 077 coupling which is model 
dependent and might be small in some models. 
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II. EVOLUTION OF THE AXION DENSITY 



In the cooling of the universe, as soon as the temperature T ~ is reached, massless 
axions can be produced by means of reactions 1), 2), and 3). The production rate has to 
compete with the expansion of the universe, characterized by the Hubble expansion rate. 



We will assume that, in the range of temperatures where we will apply the evolution 
equation for interating axions, the particle content corresponds to SU (3)c x SU (2) ^ x ?7(l)y, 
with 3 families and one Higgs boson, with all particles at the same temperature. Then, the 
relativistic degrees of freedom appearing in (|^) are given by 

7 

9* ~ Qhoson ~\~ o dfermion — 106.75 (8) 
O 

The behaviour of the axion density Ua is given in terms of the Boltzmann equation 
i 1,0, 



n T) 

-^ + 3Hna = T[n:'^-na] (9) 
Here, F is the thermal averaged interaction rate for the process a + i ^ 1 + 2 

Sv2E^2E,^ j ^^^^(2^)U*(P1+P2-P.-P.)J2\^1\' (11) 

with dp = d^p/{2'K)^. The phase space occupancy in kinetic equilibrium /"^"^ is given by the 
Fermi-Dirac or Bose-Einstein distributions, 

where we will take the relativistic limit where the chemical potential /x goes to zero. In eq 
(p!T|), |M|^ is summed over all the degrees of freedom. The equilibrium densities appearing 
in ([T0|) are given by 

n'' = gldpr (13) 



where g is the degrees of freedom of the particle. In the period of interest, thermal axions 
are massless and have a density 

n- = ^T^ (14) 



where ({3) = 1.20206 ... is the Riemann zeta function of 3. 



We follow the convenient procedure of normalizing particle densities to the 

entropy density, 



so that we define 



2-71^ 

s = —g,sT^ = OAAxg,sT' (15) 

r = ^ (16) 

During the epoch in which we will apply the Boltzmann equation (P), we can use g^s = g* = 
constant, and given by (§). Then we have that 

ye,^!C^ 45C(3) 1^0^ 
s 27r4 g^ g^ 

is constant. We can write the Boltzmann equation (P) as 



where we have defined 



^ - Y (19) 



We will now show that equation ([18|) can be integrated in our case. First, we advance 



that r oc (see sec. PD so that T/H oc 1/x and we can define the constant 



k^x^ (20) 



Second, we use the fact that Y'^'^ is independent of x. We define 



_ Y 



and, finally, we write equation ([18|) in the form 



x^^ = k{l-r]) (22) 
that has the solution 

7]{x) = 1 - e^(i/"-^) (23) 

Since we start at a; = 1 (T = F^,) with no axions, in (|23|) we have specified the initial 
condition ri{x = 1) = 0. 

From the solution (^) we see that i] starts to grow at x = 1. At x = k (that is, F = H) 
axions decouple from the QCD plasma. We call Y^ the value of Y at decoupling. For x > k 



the value of t] remains constant, so we have Y = Yd. We are interested in the situation that 
axions have in practice a thermal spectrum, so we will ask that Yd differs from Y'^'^ in less 
than 5%, 

2^ = = A;) = 1 - e^^^/''-^^ > 0.95 (24) 

which implies 

k>4: (25) 

When k satisfies this inequality, a thermal population of axions is born in the early universe. 
In the next section we calculate F, which, as we see from (^), will give us the clue of which 
is the corresponding range of Fa- 



III. CALCULATION OF THERMALLY AVERAGED INTERACTION RATES 

As expression (0) is Lorentz invariant, we can evaluate it in any reference system. For 
convenience, we choose the center of mass system and will consider that all the particles are 
massless because we are at very high energies. So, we will write (0) as 

where dcM is the usual total cross section in the center of mass system, but with no average 
over the initial degrees of freedom, and s is the Mandelstam invariant. From the Lagrangian 
piece (|^) one can easily find the Feynman rules for the couplings between axions and gluons 
that appear in the diagrams in Figs.|l|, and ^, which are necessary to calculate the cross 
sections for the considered processes. Some of these cross sections diverge logarithmically 
in the t and u channel. The processes we consider take place in a plasma with vanishing 
global color and so color is effectively screened for distances bigger than m^^, with niD the 
QCD Debye mass given by [Tl[ 



mD=9s ^ — T =8nasT (27) 

where iVc = 3 is the number of colors and Nf = 6 is the number of flavors. We cutt-off the 
divergences using the Debye mass. Cross sections have the form 

acM = A\n(^]+B (28) 



m 



D 



where A and B are constants, which for the considered processes have the numerical value 



1) ^=0 , 5=ai 
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2) A = ^|t , 5 



q^l /I — J5_ ^ R — __55. 

^ — 27r2 F? ' ^ ~ 87r2 F? 



Putting (^) in (p6D we find integrals of tlie type 

1 /•°° 

W) Jo '^"'e-±l 



^ dx-^ (29) 



, 1 /•°° , Inx 

^W^J dx-—T (30) 



C(3) 7o ± 1 

tliat can be done analitically. We liave 

I- = I (31) 

J- = 2 (32) 

L+ = i(9 + ln4-67 + 6^) (33) 

^- = 3-27 + 2^ (34) 

where 7 = 0.5772 ... is tlie Euler's Gamma constant and C'(3) = —0.1981 . . .. Witli tliese 
results, the integral (^6]) gives 

r = {[5-^(1/2 + ln(27ra,) )] I^r + AL^I~ + A J± L" } (35) 

where the sign ± depends on whether the particle that goes whith the axion is a gluon (— ) 
or a quark (+). Thus, for each of the processes 1), 2), and 3) we get, respectively, 

^ ria + g^q + q) = ^^T'Nf^-^ (36) 
12 = T{a + q^g + q) + T{a + q^g + q) = 

4 2Nf ^(L^ + h--^-^n (2vra.) - ^-^) (37) 

T, ^ r{a + g^g + g) = ^^T'^^ - In (2vra.) - (38) 

r is the sum of all Fj, and turns out to be 

r~ 7.1x10-6— (39) 

a 

where we have introduced as — 1/35, corresponding to energies E ~ 10^^ GeV. 
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IV. CONCLUSIONS 



Using (|39|), d^) and the definition ( ]20|) we have 



Fa r 15 GeV , , 

k=^-c:^5.0x 10'' -y- (40) 



a 



The inequahty ( pSD translates into 

F„ < 1.2 X 10^2 GeV (41) 

This is our main resuh. Models with Fa satisfying (|41| ) have axions thermalizing in the early 
universe and predict that today there is a thermal population of axions. 

We get a much higher value than the one in (|^). One of the reasons is that the interaction 
rate (^) increases as at high energies, while the rate coming from the photoproduction 
process considered in ^ goes as T since the axion is attached to the fermion line with 
the usual derivative coupling. The other reason is that color and flavor factors make the 
interaction rate larger. 

Our result in ( PD is in fact conservative. This is due to the fact that the effective theory 
below the scale Fa has couplings of the axion to the electroweak gauge bosons W^, W^, 
and B, similar to (|l|). This generates processes like 1), 2), and 3) with gluons replaced by 
appropiate combinations of W^, W^, and B. These new processes contain in general quarks 
and also leptons. The form of the couplings is fixed but their magnitude is model dependent 
since it is different for different PQ-charge assignments to the matter fields. Given this 
model dependence, we do not include them in our analysis. They would lead to an increase 
of the numerical value in (^9[) by a factor of ~ 2 and consequently would make the range of 
Fa in (^) larger by the same amount. 

Once we assume that axions thermalize, it is interesting to know the density of this 
thermal axion population. After decoupling, axions redshift freely until today. Their density 
today, Uao, can be calculated using the fact that after decoupling, the value of Y stays 
constant. We need the entropy density today, sq, which has the contribution of photons, 
neutrinos and axions. Even in the case that axions are relativistic today, the contribution to 
So can be neglected since the axion temperature is much colder than the cosmic microwave 
background temperature, Tq = 2.75 K. Indeed, since the time of axion decoupling, the 
decoupling of other particles have heated the photon bath but not the axion one. One gets 



that the entropy density today corresponds to g^s = 43/11 = 3.91 in (p^) 



Then, we use (H), (|T5D, and (0), to calculate the axion number density today 

27 

riao = - (0.44 X 3.91 T^) ~ 7.5 cm-^ (42) 
lUo.75 

This small density was also found by Turner ||^. At the view of the small mass (^ one can 
conclude that it is very unlikely that it can play a direct role in the efforts to detect axions 

i- 

It is also interesting to determine the range of temperatures where axions interact with the 
QCD plasma. For a given Fa satisfying (0), it is not difficult to see that the thermalization 
range is 

Fn > T > f43l 

~ 5 X 1012 GeV ^ ' 

For example, for the highest possible value of Fa where still a thermal population arises. 

Fa = 1.2 X 10^2 GeV, we get a temperature range 

1.2 X 10^2 GeV > T > 0.3 X 10^^ GeV (44) 

For the lowest value of Fa allowed by observation (see eq(H)), the range is 

6 X 10® GeV > T > 7 X 10^ GeV (45) 

Our final comment refers to a related effect also happening in the thermalization period. 
If there are non-thermal axions produced at this period, they would thermalize since they 
would interact with the QCD plasma through reactions 1), 2) and 3). Thus, no matter which 
was the energy spectrum and how high was the production rate of the originally non-thermal 
axions, they would end up with the thermal spectrum ( p!2[ ) and density (0). There is no 
effect on axions from the misalignment mechanism and domain wall decay, since they are 
generated at T ~ Aqcd, that cannot be in the range (^). Strings start to decay into axions 
after their formation at T ^ Fq. Current models have the bulk of axions produced also 
at T ~ Aqcd- If this is the case, then thermalization effects on string axions are small. 
However, the physics of axion strings is still not completely settled and we would like to 
stress that if a string model or any other source of non-thermal axions have substantial axion 
production in the range (|4^) then those axions would thermalize and finally we would have 
a single thermal population that dilute away with the expansion of the universe as we have 
studied in this paper. 
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FIG. 1: Feynman diagram for the process a + g q + q. 




FIG. 2: Feynman diagrams for the processes a + q ^ g + q and a + q'^ g + q. 




FIG. 3: Feynman diagrams for the process a + g ^ g + g. 
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